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PREFACE 


This  report  summarizes  the  work  completed  during  the  second  year  of  a  project  focused 
on  improving  the  WBMOD  model  of  ionospheric  scintillation.  This  'work  is  part  of  a  larger 
effort  which  has  the  overzdl  objective  of  providing  the  USAF  Air  Weather  Service  with  an  im¬ 
proved  capability  for  providing  support  to  their  customers  who  are  impacted  by  ionospheric 
scintillation  effects. 

We  would  like  to  theink  Dr.  Smtimay  Basu  of  the  USAF  Phillips  Laboratory  (PL/GPIA) 
for  his  guidance  and  suggestions  for  improving  the  model,  and  Dr-  Fred  Rich  (PL/ GPS  G) 
for  providing  the  DMSP  SSJ  particle  boundary  data. 


vli 


1.  Introduction 


Many  modern  military  systems  used  for  communications,  command  and  control,  naviga¬ 
tion,  and  surveillzmce  depend  on  reliable  and  relatively  noise-free  transmission  of  radiowave 
signals  through  the  earth’s  ionosphere.  Small-scale  irregularities  in  the  ionospheric  density 
can  cause  severe  distortion,  known  as  radiowave  scintillation,  of  both  the  amplitude  and 
phase  of  these  signals.  A  basic  tool  used  in  estimating  these  effects  on  systems  is  a  computer 
program,  WBMOD,  bsised  on  a  single-scatter  phase-screen  propagation  model  and  a  number 
of  empirical  models  of  the  global  morphology  of  ionospheric  density  irregularities.  A  recent 
validation  of  the  WBMOD  model  showed  that  it  was  deficient  in  a  number  of  areas.  Much 
of  this  deficiency  could  be  traced  to  the  limited  coverage  of  the  data  used  in  developing 
the  model  (a  single  station  at  high  latitudes  and  two  stations  at  equatorial  latitudes).  The, 
objective  of  this  study  is  to  investigate  and  implement  improvements  to  the  WBMOD  model 
based  on  more  extensive  data  sets  in  both  the  equatorial  and  high-latitude  regimes. 

2.  Equatorial  Model  Development 

The  development  of  the  new  equatorial  model  was  completed  during  the  report  period, 
and  the  new  model  was  implemented  in  version  12.00  of  WBMOD.  The  work  focused  on 
completion  of  the  variation  of  equatorial  scintillation  with  SSN  and  with  geomagnetic  activity 
levels. 

2.1.  Solar  Cycle  Variation.  Figure  1  shows  the  distribution  of  log(CfcZ-)  calculated 
from  VHF  observations  in  the  nighttime  (near  2200  local  time)  Ancon  and  Kwajalein 
WIDEBAND  data  sets  for  three  ranges  of  SSN,  and  Figure  2  is  a  scatter-plot  of  the  same 
data  as  a  function  of  SSN.  A  number  of  different  behaviors  can  be  noted  in  these  plots: 

1.  The  log(CjtL)  values  at  the  center  of  the  peahs  in  the  distribution  (/x  in  the  model 
PDF)  increase  with  SSN.  This  behavior  was  noted  in  earlier  modeling  [.Secan,  1989] 
and  was  included  in  previous  version  of  WBMOD. 

2.  The  width  of  the  high-5S/V  peak  (o-o.mjfcf)  decreases  with  increasing  SSN. 

3.  The  population  in  high-5SA''  peak  {rj2, night)  increases  with  SSN,  but  exhibits  what 
appears  to  be  a  saturation  effect  at  high  SSN. 


Since  the  other  data  sets  were  solely  intensity  data,  they  were  not  useful  for  determining 
the  characteristics  of  the  shape  (cr)  and  location  (fi)  of  the  PDF  peaks  due  to  distortions  at 
both  ends  of  the  distribution  (nojse-floor  and  data- reduction  induced  distortions  at  the  lower 
end  of  the  distribution,  and  saturation  effects  caused  by  Fresnel  filtering  at  the  higher  end). 
We  therfore  decided  to  use  the  WIDEBAND  ph2ise-scintillation  data  shown  in  these  figures 
to  specify  the  SSN  variation  of  these  PDF  parameters  and  to  use  the  larger  5'/(cfB)-based 
data  sets  to  specify  the  SSN  variation  of  the  parameter  which  determines  the  population 
in  the  “plume”  peak,  r)2,night-  The  cr2,night  parameter  varies  linearly  with  SSN  from  1.0  at 
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Ancon  and  Kwajalein  -  WIDEBAND  (VHF) 


26.0  28.0  30.0  32.0  34.0  36.0 

logCCJdL)  [derived  from  PRMS] 


Figiire  1:  Distribution  functions  of  log(Cibij)  calculated  from  the  Ancon  amd  Kwajalein  data 
sets.  The  log(Cjk£)  values  were  calculated  from  the  VHF  observations,  and  only  data 
from  the  night  passes  were  included.  Three  distributions  are  plotted,  one  for  each  of  three 
SSN  ranges:  SSN  <  50,  50  <  SSN  <  100,  and  SSN  >  100. 

Ancon  and  Kwajalein  -  WIDEBAND  (VHF) 


SSN 

Figure  2;  Scatter  plot  of  log(CfcL)  calculated  from  the  Ancon  and  Kwajalein  data  sets.  The 
log(C*Z)  values  were  calculated  from  the  VHF  observations,  cind  only  data  from  the 
night  passes  were  included. 
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SSN  =  0  to  0.8  at  SSN  =  SSN^,  where  SSN,  is  the  point  at  which  the  SSN  variation 
“saturates”  (this  value  is  discussed  later).  The  values  of  the  four  /i  par<iineters  also  vary 
linearly  with  SSN.  The  constants  of  the  linear  equations,  listed  in  Table  1,  were  extracted 
from  the  (r^-derived  log(C'jfel)  data  shown  in  Figure  2  and  the  corresponding  daytime  (near 
1000  local  time)  data.  Since  there  are  no  data  of  this  type  available  for  the  anomaly  crest 
region,  the  variations  of  (J^iniaht  zmd  three  of  the  four  fi  parameters  in  the  anomaly  crest 
region  are  taken  to  be  identic^  to  that  found  in  the  WIDEBAND  equatorial  data.  The  only 
exception  is  the  “DC”  term  in  the  variation  of  the  log(C*i^)  value  for  the  plume  peak  — 
this  was  adjusted  to  agree  with  the  levels  of  scintillation  observed  in  the  Ascension  Island 
database  eis  reported  in  Secan  et  ai,  [1993a]. 


Equatorial 

Anomaly  Crest 

Parameter 

Intercept 

Slope 

Intercept 

Slope 

f^\,day 

27.0 

0.0022 

27.0 

0.0022 

MS.dai/ 

0.0 

0.0 

0.0 

0.0 

29.7 

0.0002 

29.7 

0.0002 

t^i,nigkt 

33.3 

0.0022 

35.8 

0.0022 

Table  1:  SSN  variation  constants  for  the  equatorial  and  anomaly-crest  fi  models. 


The  variation  of  the  population  in  the  plume  peak  (^2)  of  the  PDF  with  SSN  was  modeled 
by  matching  the  patterns  of  the  diurnal /seasonal  variation  of  the  occurrence  statistics  for 
54  >  0.5  as  seen  in  the  SJ{dB)  data  sets.  Appendix  A  contains  the  plots  used  in  the 
modeling,  and  Table  2  summarizes  the  number  of  data  points  in  each  plot,  the  number  of 
points  for  which  S4  exceeded  0.5,  and  the  percent  of  time  that  S4  exceeded  0.5.  Figure  3 
illustrates  the  variation  of  the  54  >  0.5  occurrence  statistic  as  a  function  of  SSN.  The  data 
used  in  generating  this  figure  were  limited  to  low  geomagnetic  activity  levels  {Kpss  ^  2o, 
where  KpSs  is  the  A'p value  at  local  sunset),  the  first  four  hours  after  tg  =  0,  eind  to  the  two 
sixty-day  periods  covering  the  periods  of  maximum  scintillation  activity  for  each  station. 
[Note:  Many  of  the  plots  in  this  report  are  of  the  occurrence  statistic  defined  as  the  percent 
of  the  time  that  S4  exceeds  0.5.  This  statistic  will  be  denoted  Pr54(0.5)  in  the  remainder  of 
this  report.] 

In  studying  these  plots,  we  have  found  the  following  behaviors: 


1.  Frs4(0.5)  increases  with  increasing  SSN,  with  a  tendency  to  reach  “saturation”  at  a 
SSN  value  of  125  to  150. 

2.  The  saturation  levels  for  the  two  equatorial-VHF  data  sets  are  in  fairly  good  agreement, 
but  the  level  at  Huancayo  is  significantly  higher  at  low  SSN  than  at  Manila. 
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SSN 

Ground  Station 

Range 

Ascension 

Manila 

Huancayo 

34,262 

MSSSM 

56,746 

All 

2,544 

9,995 

7% 

■ga 

18% 

10,092 

<  50 

476 

1,949 

>1% 

5% 

10,331 

12,651 

50  -  100 

1,068 

1,825 

8% 

14% 

8,615 

7,422 

13,185 

100  -  150 

1,098 

1,643 

3,171 

13% 

22% 

24% 

4,569 

2,244 

150  -  200 

706 

609 

15% 

27% 

26% 

1,709 

>  200 

274 

16% 

37% 

26% 

L-band 

VHF 

VHF 

Table  2:  Number  of  data  points  as  a  function  of  SSN.  The  three  numbers  in  ezich  cell  are 
the  totzd  number  of  o’  servations  (top),  the  niunber  of  observations  of  S4  >  0.5  (middle), 
and  the  percent  of  time  that  S4  exceeded  0.5  (bottom). 


3.  The  ratio  between  the  ne.'’T-equinoctal  maxima  and  the  local-summer  minimum  reduces 
with  increasing  SSN. 

4.  The  decay  with  time  after  the  diurnal  maximum  decreases  with  increasing  SSN  at  the 
equatorial  stations,  although  not  at  Ascension  Island. 


These  variations  have  been  included  in  the  model  via  the  PDF  parameter  that  describes  the 
nighttime  plume  population,  T}2,nigKti  ®uid  the  diurnal  variation  of  the  r]2  parameter. 

The  value  for  T]2,nigkt  is  calculated  from 


Vi, night  =  fs{SSN,4>r,)fk{Kj,ss)fd{4>Urm)  (1) 

where  <j>^  is  a  longitude  parameter,  Kpss  is  the  Kp  valid  at  local  sunset,  and  (purm  is  the 
angle  between  the  sunset  terminator  and  the  local  geomagnetic  meridian  calculated  at  the 


4 


Variation  of  %S4>0j5  with  SSN 


Figure  3:  Variation  of  Prs4(0.5)  as  a  function  of  SSN  at  Ascension  Island  (L-band),  Manila 
(VHF),  and  Huancayo  (VHF).  The  data  used  were  limited  to  low  geomagnetic  activity, 
post-sunset,  and  near-equinox  (sec  text  for  details). 

IPP.  The  three  functions  describe  the  variation  of  this  parameter  with  SSN,  Kpss  (described 
in  the  next  section),  surd  season.  The  “zero”  level  is  set  in  the  f,  function,  and  the  other 
two  functions  modulate  this  level. 

The  SSN  variation  shown  by  the  data  was  a  steady  increase  in  Tj2,night  with  SSN,  from 
a  minimum  value  at  SSN  =  0  to  a  maximum  at  some  “saturation”  value  of  SSN,  denoted 
SSNf  The  model  transitions  from  the  minimum  value  to  the  saturation  value  using  a  simple 
gaussian  function  of  SSN  as  follows: 


/*  =  + 


SSN'  -  SSN,' 


nr 

rv., 


21 


J 


(2) 


where  tjo  and  rj,  are  the  minimum  and  saturation  values  for  rf2,night,  SSN'  is  the  minimum  of 
the  observed  SSN  and  the  value  of  SSN„  and  W„n  is  calculated  from  other  parameters  so 
that  ft  becomes  arbitrarily  close  to  the  value  of  rjo  for  SSN'  =  0. 

In  fitting  this  modeled  variation  to  the  data,  it  became  apparent  that  the  SSN  variation 
at  Huemcayo  and  Manila  were  impossible  to  reconcile  without  invoking  a  longitudinal  differ¬ 
ence  between  the  two  stations.  This  was  included  in  the  model  by  applying  a  longitudinal 
modula  ion  on  the  t/o  parameter  such  that  the  value  of  this  parameter  in  the  vicinity  of 
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Manila  is  a  factor  of  0.4  times  the  value  in  the  Huancayo  sector.  The  very  non-uniform 
longitudinal  coverage  of  the  data  set  made  it  impossible  to  define  this  variation  rigorously, 
so  it  has  been  implemented  in  the  following  simple  form: 


^  =  [y}ton  +  {i-mon)fi(<i>)]rioo  (3) 

where  r}ian  is  0.4  and  /j  is  given  by 

/,W)  =  0.5  [erf  -  erf  (^)]  (4) 

where  <l>  is  the  geographic  longitude;  4>i  and  02  are  60“  and  210“,  respectively  and  A0  is  set 
to  15°.  These  values  were  set  in  order  to  move  the  treursitions  away  from  the  locations  where 
we  had  data,  and  were  based  on  review  of  the  results  documented  in  Basu  et  al.  [1976]  and 
Maruyama  and  Matuura  [1984]  from  in-situ  observations.  [Note:  This  longitudinal  variation 
replaces  the  one  described  in  Scientific  Report  #  1  for  this  project  {Secan  et  al.  [1993a], 
denoted  SR-1)]. 

The  values  for  the  constants  ijg  (in  Equation  [2])  and  rjoo  (in  Equation  [3])  were  extracted 
from  the  data  sets,  and  are  given  in  Table  3. 


Equatorial 

Anomaly  Crest 

Voo 

0.70 

0.05 

1.00 

0.80 

Table  3:  Values  for  T700  and  r/,  in  the  equatorizd  and  anomaly-crest  models. 

The  models  for  the  SSN  component  of  the  diurnal  and  seasonal  variations  are  discussed 
in  Sections  2.3  and  2.4,  respectively. 

2.2  Geomagnetic  Activity  Variation.  As  described  1  SR-1,  we  have  decided  to  use  the 
value  of  Kp  at  local  sunset  (denoted  Kpss)  for  modeling  the  variations  in  scintillation  levels 
in  the  post-sunset  period  and  the  Kp  value  at  the  observation  time  for  variations  in  the  levels 
after  the  post-sunset  peak  in  scintillation  activity.  Contour  plots  of  the  occurrence  statistic 
Prs4{0.5}  similar  to  those  generated  for  the  95/V  variation  are  provided  in  Appendix  B  as  a 
function  of  Kp  and  in  Appendix  C  as  a  function  of  Kpss-  The  number  of  data  points  used 
lo  generate  each  figure,  the  number  of  points  in  each  plot  where  S4  >  0.5,  and  the  percent 
of  time  that  S4  exceeded  0.6  are  listed  in  Tables  4  (zis  a  function  of  Kp)  and  5  (zis  a  function 
of  Kpss). 

The  steady  decreztse  in  the  occurrence  statistic  evident  in  Table  5  car.  be  seen  clearly 
in  Figure  4,  which  shov's  the  variation  of  Pr 54 (0.5)  as  a  function  of  KpSs  from  four  of  the 
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Kp 

Ground  Station 

Range 

Ascension 

Manila 

Huancayo 

75,613 

80,155 

124,474 

All 

4,437 

8,210 

18,654 

6% 

10% 

15% 

33,929 

28,621 

57,263 

<  2o 

2,407 

3,335 

8,395 

7% 

12% 

15% 

32,437 

38,145 

2+  -  40 

3,758 

7,699 

5% 

15% 

14,731 

>  4+ 

4% 

17% 

L-band 

VHP 

VHP 

Table  4:  Number  of  data  points  as  a  function  of  Kp.  The  three  numbers  in  each  cell  are  the 
total  number  of  observations  (top),  the  number  of  observations  of  S4  >  0.5  (middle),  and 
the  percent  of  time  that  54  exceeded  0.5  (bottom). 


Kpss 

Ground  Station  { 

Range 

Ascension 

Manila 

Huancayo 

75,613 

80,155 

124,474 

All 

4,437 

8,210 

18,654 

6% 

15% 

:  4,262 

31,181 

56,746 

<  2o 

2,544 

9,995 

7% 

13% 

18% 

31,743 

37,415 

51,679 

2+  -  40 

1,657 

3,428 

6,813 

5% 

11% 

13% 

9,608 

11,559 

>  4+ 

236 

579 

2% 

5% 

l-band 

VHP 

VHP 

Table  5  Number  of  data  points  as  a  function  of  Kpss-  The  three  numbers  in  ead’  cell  are 
the  tot  ,1  number  of  observations  (top),  the  number  of  observations  of  S4  >  0.5  (.middle), 
ajid  the  percent  of  time  that  S4  ex  ceded  0.5  (bottom). 


Variation  of  %  S4  >  0.5  with  Kp(SS) 


Figure  4:  The  vzu-iation  of  Pri’4(0.5)  as  a  function  of  Kpss  at  Ascension  Island  (L-band), 
Manila  (VHF),  and  Huancayo  (VHF).  The  data  used  were  limited  to  post-sunset  and  near- 
equinox  (see  text  for  details). 

data  sets.  The  data  used  in  generating  this  plot  were  limited  to  the  first  four  houns  after 
fe  =  0,  and  to  the  two  sixty-day  periods  covering  the  periods  of  maximum  scintillation 
activity  for  each  station.  With  the  exception  of  Manila,  the  data  sets  show  a  steawiy  decrease 
in  Pr54(0.5)  with  Kpssi  with  a  tendency  for  the  decre2ise  to  “level  out”  at  a  KpSs  value 
between  5°  and  6^.  A  irons  [1985]  included  a  similar  decrease  with  increatsing  Kp  in  his 
model  for  the  variation  of  SI{dB)  at  Huancayo.  Figure  5  shows  the  variation,  of  SI{dB)  as 
a  function  of  time  for  four  values  of  Kp. 

We  have  modeled  this  variation  by  varying  the  Ti2,night  parameter  as  a  function  of  Kpss 
as  follows: 


n--^1.0-0.UK;gs  (5) 

where  A^'55  is  the  minimum  of  the  observed  Kpss  and  a  value  of  5.0.  The  value  of  0.14 
in  Equation  [5]  fit  both  the  Hucincayo  and  Ascension  Island  data  well,  although  the  Manila 
data  showed  a  slower  drop  with  Kpss  for  low  activity  levels. 

In  studying  the  various  plots  of  changes  in  the  diumal/seasonal  v^iation  of  Prs4(0.5) 
with  geomagnetic  activity  we  have  found  the  following: 

1.  The  ratio  of  the  scintillation  level  at  the  seasonal  maxima  (as  measured  by  Pr54(0.5))  to 
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Huancayo  (VHF)  [Aarons  Model] 


Figure  5:  Variation  of  Sl{dB)  as  a  function  of  time  for  four  levels  of  geomagnetic  activity 
taken  from  the  model  of  Aarons  [1985]. 

the  level  at  the  local  summer  minimum  tends  to  decrease  with  increasing  geomagnetic 
activity  at  the  equatorial  stations,  but  not  at  Ascension.  This  behavior  was  also 
included  in  the  Aarons  model,  as  can  be  seen  in  Figure  6. 

2.  As  the  level  of  geomagnetic  activity  increases,  the  value  of  Fr54(0.5)  increases  in  the 
late  evening/eariy  morning  hours.  This  car  cilso  be  seen  in  the  Aarons  model  (see 
Figure  5). 


These  behaviors  have  been  included  in  the  diurnal  and  seasonal  variation  models  described 
in  the  next  two  sections. 

2.3  Diurnal  Variation.  The  model  of  the  diurnal  variation  described  in  SR~1  ha-s  been 
modified  to  include  SSN  and  geomagnetic  effects  described  in  the  previous  two  sections. 
The  parameter  that  specifies  the  time  after  sunset  at  which  ‘■he  peak  level  of  scintillation  is 
reached,  t,,,  has  been  modified  to  vary  with  SSN  as  follows; 


tij  —  iatO 


fS^'\ 

\SSNs) 


(6) 


where  SSN'  and  S5/V,  are  as  defined  in  Section  2.1.  The  p<irameter  which  controls  the  rate 
at  which  the  scintillation  level  decreases  after  the  post-sunset  maximum,  W^r,  now  varies 


9 


Huancayo  (VHF)  [Aarons  Model] 


Day  of  the  Year 

Figure  6:  Variation  of  Sl{dB)  as  a  function  of  day  of  the  year  for  four  levels  of  geomagnetic 
activity  takeu  from  the  model  of  Aarons  [1985]. 

with  both  SSN  and  Kp  as  follows: 


where  is  the  minimum  of  the  observed  Kp  and  a  value  of  5.0. 

In  a  separate  study,  we  found  discontinuities  in  the  diurnal  variation  of  log(  C*Z-)  in  the 
new  equatorial-regiou  models.  We  have  modified  the  way  in  which  the  diurnal  variation  is 
implemented  in  order  to  remove  these  discontinuities.  These  changes  were  as  follows: 

1 .  T  he  function  used  to  model  the  diurnal  variation  for  the  tr,-  and  7]i  parameters  (Equation 


[6]  in  SR-1)  was  changed  to 


fiQ  = 


exp 

V 

«) 

1.0 

if  te<te<  tsr 

exp 

[-  (He)’] 

if  <  te  <  tc 

exp 

[  i‘^)  1 

exp  - 

if  te  >  tc 

(8) 


where  is  the  time  past  sunset,  t,,  and  W,,  sure  the  peak  time  and  half-width  ol 
the  post-sunset  gaussian,  t„  and  are  the  time  at  which  the  late-evening/early- 
moruing  “sunrise”  segment  of  the  variation  starts  its  decrease  and  the  half-width  in 
this  segment,  and  tc  and  Wc  are  the  time  at  which  the  late-morning  segment  of  the 
variation  starts  and  the  half-width  in  this  segment. 

2.  The  half-width  is  calculated  from 

W,,  =  0.466<„  (9) 

rather  than  a  fixed  value  of  0.80'‘. 

3.  The  t,r  time  is  calculated  from 

t,,  =  t„-f- 1.5^  (10) 

rather  than  a  fixed  value  of  1.75*. 

4.  The  Wc  half-width  in  Equation  [8]  is  calculated  from 


iyc  =  0.466  (to --tc)  .  (11) 

•5.  The  diurnal  variation  for  the  p,  parameter  has  been  altered  so  that  the  (non-plume 
population)  and  ^2  (plume  population)  parameters  have  slightly  different  variations. 
The  new  variations,  which  replace  Equation  [8]  from  SR-1,  are 


f  fUc)  if  U  <  t,r 

/^n-p(tc)  =  \  if  4.  <  <«  <  <0 

0  if  te  >  to 


(12) 
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for  the  non-plume  population,  and 


1.0 

U(Q  =  (13) 

.  /(ie)  if  t,  >  tjr 

for  the  plume  population. 

These  modifications  do  not  change  measurably  the  results  presented  in  SR  -1,  as  the  dis¬ 
continuities  only  appear  for  low  scintillation  levels  (^4  <  0.1)  and  low  frequencies  (<  100 
MHz). 

The  model  constants  in  Equations  [6]  and  (13j  arc/given  in  Table  6. 


Table  6:  Values  for  the  constants  in  the  diurnal  variation  for  the  equatorial  and  anomaly-crest 
models. 

2.4  Seasonal  Variation.  The  model  for  seasonal  variation  has  been  somewhat  changed 
from  the  description  in  SR-1.  The  form  of  the  variation  has  not  changed,  but  the  way  in 
which  the  W,  parameter  is  calculated  has.  Rather  than  “summer”  and  “winter”  values,  the 
modeling  for  W,  is  in  terms  of  “high”  and  “low”  values.  These  correspoud  to  the  values  to 
be  used  in  the  solstice  where  the  scintillation  levels  are  either  high  (with  respect  to  the  levels 
at  the  other  solstice)  or  low.  The  determination  of  which  value  tor  W,  is  valid  is  a  function 
of  the  season  and  of  the  lonj^itude  sector. 

The  “low”  value  for  li',,  Wsjowi  corresponds  to  what  had  been  called  the  “winter”  value, 
which  has  been  changed  to  a  fixed  value  of  13.0.  The  “high”  value,  corresponding 

to  the  old  “summer”  value,  now  varies  with  both  SSN  and  Kpss  in  order  to  replicate  the 
seasonal  behavior  seen  in  the  data.  This  parameter  is  now  calculated  from  the  follorving: 


where  the  values  for  the  three  constants  in  this  equation  for  the  equatorial  eind  amomaly-crest 
models  are  given  in  Table  7. 


Equatorial 

Anomaly  Crest 

16.0 

16.0 

-0.01 

-0.01 

1.0 

0.0 

Table  7:  Values  for  the  constants  in  the  W,^high  equation  for  the  equatorial  and  anomaly-crest 
models. 

The  selection  of  which  of  the  two  values  to  use  for  W,  is  based  on  the  sezison  and  on  the 
longitude  sector.  In  the  Pacific  sector  (where  //  =  1.0  (see  Equation  [4])),  Wg^higk  is  used 
on  the  June  solstice  side  of  the  seasonal  activity  peaks  and  is  used  on  the  December 

solstice  side.  In  the  Atlantic  sector  (//  =  0.0),  is  used  on  the  June  solstice  side 

of  the  seasonal  activity  peaks  and  W^high  is  used  on  the  December  solstice  side.  This  is 
implemented  in  the  model  code  as  follows: 


(15) 


where 


m  = 


< 


//(» 

if  <  0 

i-m) 

if  ^term  ^  0 

(16) 


The  rationale  for  this  aspect  of  the  model  is  the  obvious  difference  in  the  solstitial  behavior 
at  the  three  Atlantic  sector  stations  (Huancayo,  Ancon,  and  Ascension  Island)  and  that  seen 
at  the  two  Pacific  sector  stations  (Manila  and  Kwajalein).  This  difference  has  been  noted 
before  (see  Basv.  et  al.,  [1976]  and  Maruyama  and  Matuura,  [1984]),  and  was  included  in  a 
soon  vhat  crude  fetshion  in  the  last  upgrade  of  the  WBMOD  equatorial  model  [Secan,  1989]. 

2.5  WBMOD  Version  12.  The  new  models  for  equatorial  scintillation  described  in  this 
report  and  SR-1  have  been  implemented  in  Version  12  of  WBMOD.  This  version  of  WBMOD 
has  the  following  changes  from  Version  1111: 


1.  The  model  for  the  irregularity  strength  parameter,  CkL,  has  been  expanded  from 
modeling  merely  the  mean  CkL  value  to  modeling  the  full  log(CitL)  PDF. 

2.  The  equatorial  model  for  CkL  has  been  replaced. 


1.3 


3.  The  cailculation  of  S4  on  a  two-way  path  bztsed  on  the  theory  developed  by  Fremouw  and 
Ishimaru  [1992]  has;  been  implemented.  The  user  is  requested  to  enter  the  correlation 
between  the  up-  and  down-links,  a  parameter  which  will  be  calculated  by  the  model  in 
a  later  implementation. 

4.  The  user  is  given  two  options  for  model  output:  (1)  the  occurrence  statistics  for  a 
user-specified  level  of  scintillation,  or  (2)  the  level  of  scintillation  for  a  user-specified 
occurrence  percentile.  The  user  may  specify  the  scintillation  level  either  in  terms  of  S4 
or  of  <7^  (if  17^  is  to  be  calculated). 

5.  The  SI{dB)  intensity  scintillation  parameter  is  now  calculated  and  included  in  the 
model  output. 

A  User  Guide  wais  written  for  the  reseawch  version  of  WBMOD.  It  presents  a  brief  de¬ 
scription  of  the  model  and  describes  how  the  model  should  be  run  [5ecan,  1993].  Included 
in  the  User  Guide  are  copies  of  papers  presented  at  the  1987  and  1993  Ionospheric  Effects 
Symposia  [5ccan  et  al.,  1987;  Secan  et  ai,  1993b]  tc  augment  the  model  description.  Copies 
of  this  document  were  provided  to  PL/GPIA  and  to  the  Air  Force  Space  Forecast  Cen¬ 
ter  (AFSFC).  The  AFSFC  also  received  a  copy  of  the  reserirch  version  of  WBMOD  12.02 
compiled  to  run  on  an  Intel  80386-class  computer. 

3.  High’Latitude  Model  Development 

With  completion  of  the  equatorial  region  of  the  model,  attention  has  been  moved  to 
upgrading  the  high-latitude  region.  A  very  different  set  of  data  will  be  used  in  this  part  of 
the  model-upgrade  effort  from  that  used  in  the  equatorial  upgrade.  As  with  the  equatorial 
database,  one  of  the  components  of  the  high-latitude  database  is  a  set  of  data  taken  from  the 
WIDEBAND  satellite  experiment,  collected  at  Poker  Flat,  Alaska.  The  remaining  data  were 
collected  from  the  HiLat  and  Polar  BEAR  satellites  from  four  ground  stations:  Bellevue, 
Washington;  Tromso,  Norway;  Ft.  Churchill,  Canada;  and  Sondre  Stromfjord,  Greenland. 
Table  8  lists  information  about  each  of  these  five  data  sets.  The  column  labeled  IPP  in  this 
table  is  the  geomagnetic  (apex)  latitude  at  the  350-km  ionospheric  penetration  point,  and 
is  the  average  for  each  data  set. 

3.1  Model  Plan.  I'heie  will  be  six  main  stages  to  the  high-latitude  upgrade  effort: 


1 .  Establish  the  software  and  procedures  for  calculating  estimates  of  CkL  from  the  various 
scintillation  observations.  This  has  been  completed  at  all  stations,  but  we  are  concerned 
about  the  potential  effects  of  small  errors  in  the  satellite  locations  in  the  HiLat  and 
Polar  BEAR  modeling  databases  on  the  calculated  CkL.  These  errors  are  due  to  the 
use  of  outdated  ephemeris  sets  at  field  sites  at  times  when  updated  ephemerides  were 
not  available.  This  is  discussed  further  in  Section  3.2. 

2.  Update  the  geophysical  variations,  i.e.,  the  variation  of  CkL  with  SSN  and  Ky. 
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Ground  Station 

Satellite 

IPP 

Dates 

Points 

Bellevue 

Poker  Flat 

Tromso 

Ft.  Churchill 

Sondre  Strom^ord 

HiLat/PB 

WIDEBAND 

HiLat/PB 

HiLat/PB 

HiLat/PB 

54.1® 

65.9® 

67.1® 

69.7® 

74.0® 

1984-1937 

1976-1979 

1984-1987 

1984-1988 

1984-1988 

86,240 

61,340 

512,490 

317,090 

696,020 

Data  types  — 

Si  :  Intensity  scintillation  index  (RMS) 
:  Phase  scintillation  index  (RMS) 

T  :  Phase  PDS  spectral  strength  (iHz) 
p  :  Phase  PDS  spectral  slope 


Table  8;  Data  available  for  use  in  the  high-latitude  model  upgrade. 


3.  Update  the  diurnal  variation. 


Update  the  seatsonal/longitudinal  variation.  The  present  WBMOD  does  not  include  a 
seasonal  or  longitudinal  variation  as  none  was  found  in  the  WIDEBAND  Poker  Flat 
data  set.  There  is  evidence  of  a  definite  seasonal  variation  at  other  longitude  sectors, 
a  variation  which  has  been  tied  to  the  offset  between  the  geographic  and  geomagnetic 
poles.  We  plan  to  use  the  longitude  diversity  in  our  data  set  to  address  this  issue. 


5.  Implement  all  changes  to  th?  high-latitude  model  into  the  research  version  of  WBMOD. 
This  version  will  be  denoted  Version  13.00. 


3.2  CkL  Calculation,  Crlculating  CkL  from  the  scintillation  observations  in  the  high- 
latitude  database  is  done  using  the  algorithms  developed  for  the  Air  Weather  Service  (AWS) 
GBCKL  program  [5ecan  ef  al.,  1990],  These  are  essentially  the  algorithms  used  in  WB¬ 
MOD  to  calculate  .scuitillation  indices  from  CkL  run  backwards.  These  algorithms  require 
specification  of: 


1.  Details  of  the  ray  path  (t.c.,  the  satellite- to-ground  geometry).  This  is  calculated  using 
straightforward  spherical  trigonometry  algorithms. 


2.  Details  of  the  orientation  and  shape  of  the  plasma-density  irregularities  which  cause 
the  scintillation.  This  information  is  not  available  from  observations,  so  it  is  calculated 
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from  the  models  in  the  current  WBMOD  for  the  axial  ratios  a  zmd  b  and  for  the 
orientation  angle  S. 

3.  The  velocity  of  the  line-of-sight  with  respect  to  the  plasma-density  irregularities.  There 
are  two  components  to  this:  velocity  due  to  motion  of  either  end  of  the  ray  path  (satel¬ 
lite  motion,  for  example),  and  the  in-situ  drift  velocity  of  the  irregularities.  The  com¬ 
ponent  due  to  the  satellite  motion  is  calculated  using  further  spherical  trigonometry, 
and  the  in-situ  drift  is  modeled  using  the  drift-velocity  model  in  the  current  WBMOD. 


Estimates  of  CkL  can  be  calculated  from  either  S4  and  p,  cr^  gmd  p,  or  T  and  p.  For  the 
modeling  effort,  we  will  calculate  CkL  from  and  S4  with  p  from  the  WIDEBAND  data, 
set  auid  from  T  with  p  from  the  HiLat  and  Polar  BEAR  data  sets. 


As  stated  earlier,  there  may  be  a  problem  in  the  CkL  values  due  to  errors  in  the  satellite 
locations  stored  in  the  modeling  database  files.  This  can  cause  severe  problems  when  the 
line  of-sight  is  passing  very  nearly  parallel  to  geomagnetic  field  lines  or  nearly  tangent  to 
L-shells.  Scintillation  tends  to  be  enhanced  at  these  points,  and  small  errors  in  calculating 
the  angle  between  the  line-of-sight  and  the  direction  of  the  geomagnetic  field  can  result  in 
non-trivial  errors  in  the  calculated  CkL.  We  will  be  looking  into  this  over  the  next  period 
by  recalculating  all  satellite  locations  (HiLat  and  Polar  BEAR)  using  a  Isu^ge  database  of 
ephemerides  for  these  two  satellites  provided  by  Dr.  Kelso  at  the  US  Air  Force  Institute  of 
Technology. 

3.3  Preliminary  Results. 

3.3.1  Scintillation  Boundary.  The  initial  plan  for  modeling  the  transition  boundary  from 
the  mid-latitude  to  auroral  regions  was  to  derive  the  location  and  behavior  of  this  boundary 
directly  from  the  scintillation  observations.  This  has  been  found  to  be  impossible  to  do  using 
the  large  binned  data  sets  we  routinely  use  in  the  modeling  work.  The  boundary  becomes 
smeared  in  both  latitude  and  time  when  many  passes  are  averaged  together.  There  are 
obvious  changes  in  the  average  patterns  (plots  of  average  log(CfcL)  as  a  function  of  latitude 
and  local  time)  with  Kp.  but  they  are  very  small  and  difficult  to  model.  The  alternatives  to 
extracting  the  boundary  behavior  from  the  aggregate  database  are  (1)  to  manually  inspect  a 
large  number  of  individual  passes  to  determine  where  the  boundziry  ' ;  located  in  each  pass, 


aggregate  database  and  determine  boundary  locations  for  each  pass.  Neither  of  these  tasks 
is  within  the  scope  of  ti  e  present  work. 

As  a  workable  alternative,  we  are  going  to  follow  the  path  used  in  developing  the  original 
WBMOD  high-latitude  model  and  use  a  boundary  derived  from  another  data  source  as  the 
basis  for  the  model  of  the  scintillation  boundary.  The  boundary  of  choice  is  that  specified  by 
the  DMSP  SSJ/4  electron- precipitation  data,  and  a  model  based  on  those  data.  As  a  first 
look.,  we  looked  at  the  behavior  of  ing(CfcZ/)  as  a  function  of  the  location  of  a  precipitation 
boundary  derived  from  a  model.  This  was  promising,  so  the  full  SSJ/4  boundary  databause 
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was  obtained  from  Dr.  Fred  Hicti  (PL/GPSG)  to  include  in  our  modeling  database. 

Our  choice  of  the  SSJ  precipitation  boundary  was  based  on  several  factors.  First,  there 
is  a  large  database  avEiilable  of  this  boundary  location  covering  most  of  our  scintillation 
database.  Second,  there  s  a  model  available  for  use  in  WBMOD  which  will  provide  this 
boundary  location  as  a  function  of  K^,.  Finally,  this  is  a  parameter  which  is  routinely 
available  at  the  AFSFC. 

The  SSJ  model  used  is  based  on  offset  circles  in  geomagnetic  latitude  and  local  time 
fitted  to  the  model  of  the  boundary  movement  as  a  function  of  Kj,  presented  in  Gusstnhoven 
et  al.  [1983].  In  that  paper,  they  presented  different  linear  functions  of  the  boundary  latitude 
with  Kp  for  each  local  geomagnetic  time.  The  present  model  defines  the  boundary  location 
with  the  radius  aind  center  (latitude  and  local  time)  as  functions  of  Kp.  These  three  models 
are  as  follows: 


R.  =  20.97  +  1.658A'p  degrees 

L  =  2.466  -I-  0.3651  A'p  degrees  (17) 

=  2.862  -  0.09638AV  +  0.004893/^^*  hours 


This  model  has  been  used  in  the  WBMOD  high-latitude  drift  velocity  model  since  Version 
lOHl  to  specify  the  behavior  of  the  drift  velocity  equatorward  of  the  reversal  boundary. 

The  SSJ/4  boundary  database  provided  by  Dr.  Rich  is  the  boundary  location  interpo¬ 
lated  from  the  local  time  of  observation  to  local  geomagnetic  midnight.  This  database  covers 
the  years  1983  through  1991  (all  of  the  HiLat  and  Polar  BEAR  data  sets,  but  none  of  the 
WIDEBAND  data  set),  and  contains  over  122,000  data  points. 


To  determine  the  behavior  of  the  scintillation  boundary,  we  will  study  the  v^uiation  of 
CkL  with  respect  to  the  location  of  the  precipitation  boundary.  This  should  remove,  at 
least  to  first  order,  A'p-dependent  effects  in  the  scintillation  due  to  motion  of  the  boundary. 
[Note:  Recall  that  W'BMOD  is  a  climatological  model,  and  not  a  phenomenological  one. 
While  this  leist  statement  is  workable  for  a  climatological  model,  it  is  not  necessarily  true  for 
individual  cases.  The  scintillation  boundary  is  most  likely  less  dynamic  than  the  precipitation 
boundary.] 


We  have  just  begun  to  work  on  this  study,  but  the  preliminary  remits  are  very  promising. 
Figures  7  through  11  show  the  variation  logfC^A)  (averaged  in  one-degree  by  1-hour  bins) 
as  a  function  of  distance  from  the  precipitation  boundary  (in  geomagnetic  latitude)  and 
apex  (geomagnetic)  local  time.  The  boundary  locations  for  all  of  the  FiLat  and  Polar 
BEAR  data  were  derived  from  the  boundary  database.  Only  boundaries  observed  within 
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Figure  7;  Bin-averaged  behavior  of  log(Cfcl-)  at  Bellevue  <is  a  function  of  the  distance  in 
geomagnetic  latitude  from  the  high-latitude  precipitation  boundary  and  geomagnetic  local 
time.  The  boundary  location  was  derived  from  the  DMSP  SSJ  boundary  database. 
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Figure  8:  Bin-averaged  behavior  of  log(CjtL)  at  Poker  Flat  as  a  function  of  the  distance  in 
geomagnetic  latitude  from  the  high-latitude  precipitation  boundajy  and  geomagnetic  local 
time.  The  boundary  location  was  derived  from  the  A'p-dependent  model  described  in  the 
text. 


18 


Tromso-HL/PB(VHF) 
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X-axfs  bin  size:  1.00  Z  parameter  iog(CkL)  [firon;  S4] 

Y-axis  bin  size;  2.00  Zmin;  30.8  Z  max;  32.8 

Number  of  data  points:  181227  White:  30.5  Black:  32.5 

Data  limited  to  1.5  <  p  <  4.5 
Data  limited  to  THETA  >  65dc£ 

VDRJFT  from  WBMOD  model 

Figure  9:  Bin-averaged  behavior  of  log(Cfc£)  at  Tromso  as  a  function  of  the  distance  in 
geomaignetic  latitude  from  the  high-latitude  precipitation  boundary  and  geomagnetic  local 
time.  The  boundzwy  location  was  derived  from  the  DMSP  SSJ  boundary  database. 
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Figure  10;  Bin-averaged  behavior  of  log(CfcL)  at  Ft.  Churchill  as  a  function  of  the  distance 
in  geomagnetic  latitude  from  the  high-latitude  precipitation  boundary  and  geomagnetic  local 
time.  The  boundary  location  was  derived  from  the  DMSP  SSJ  boundary  database. 
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Figure  11;  Bin-averaged  behavior  of  log{C'jfe£)  at  Sondre  StromQord  as  a  function  if  the 
distance  in  geomagnetic  latitude  from  the  high-latitude  precipitation  boundary  and  geo¬ 
magnetic  local  time.  The  boundary  location  was  derived  from  the  DMSP  SSJ  boundary 
databzLse. 
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30  minutes  of  the  scintillation  observation  were  used.  The  boundairy  location  for  the  Poker 
Flat  (WIDEBAND)  data  set  w<is  calculated  using  Kp  and  the  SSJ  boundary  model.  All 
of  these  plots  show  systematic  behavior  with  respect  to  the  precipitation  boundary,  which 
shows  promise  for  the  modeling  effort. 

3.3.2  Kj,  Variations.  The  current  high-latitude  CkL  model  moves  the  location  (latitudinal 
extent)  and  transition-width  (“steepness”  of  the  transition)  of  the  scintillation  boundary  as 
a  function  of  Kp.  This  boundary  moves  equatorward  with  Kp  and  the  width  increases 
accordingly.  The  Ck  L  value  at  a  fixed  station  increases  with  increasing  Kp  solely  due  to  the 
motion  of  the  boundary.  The  level  of  CkL  within  the  high-latitude  region  does  not  change 
as  a  function  of  Kp. 

Our  preliminary  study  shows  that  this  is  in  rough  agreement  with  the  present  data  set, 
although  we  do  see  some  increase  in  CkL  with  increaising  Kp  for  locations  which  are  well 
within  the  boundary  location.  We  will  be  looking  closely  at  this  to  determine  if  this  is  a  real 
effect,  cr  a  residual  effect  of  errors  in  locating  the  boundary. 

3.3.3  SSN  Variations.  All  of  the  data,  sets  .show  a  gener?!  incrccise  in  scintillation  levels 
with  increasing  SSN.  This  is  not  surprising  since  the  overedl  ionospheric  density  is  known 
to  increase  with  increasing  SSN.  Our  initial  results,  however,  show  that  the  increase  of  CkL 
with  SSN  is  slower  than  that  in  the  present  WBMOD.  The  current  model  shows  log(CT£) 
increasing  as  0.017  times  the  SSN,  while  the  initial  analysis  of  the  data  at  hand  shows  only 
an  increase  of  0.005  times  the  SSN.  Re-analysis  of  the  WIDEBAND  data  from  Poker  Flat 
(which  was  the  only  data  set  used  in  the  original  modeling)  agrees  with  this  new  value  for 
the  variation  with  SSN.  We  are,  at  this  time,  at  a  loss  to  explain  this  discrepancy.  We  have 
seen  evidence  for  a  slightly  different  variation  with  SSN  in  different  latitude  and  local-time 
regimes.  We  will  investigate  this  further. 

4.  Conclusion 

The  upgrade  to  the  equatorial  region  of  the  WBMOD  CkL  model  has  been  completed, 
and  the  upgrade  has  been  implemented  in  the  research  version  of  WBMOD  (Version  12).  A 
User’s  Guide  has  been  written  which  provides  a  brief  description  of  the  model,  and  describes 
how  the  model  should  be  run  and  the  inputs  requested  by  the  model  as  it  runs. 

V,  >rk  has  begun  on  the  upgrade  to  the  high-latitude  CkL  model.  This  upgrade  will  use 
data  from  the  WIDEBAND  (Poker  Flat),  and  HiLat/Polar  BEAR  (Bellevue,  Tromso,  Ft. 
Churchill,  Sondre  Stromfjord)  experiments.  When  the  high-latitude  model  is  completed,  the 
new  algorithms  will  be  implemented  into  the  research  version  of  WBMOD  (Version  13). 
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Appendix  A,  Data  Plots  —  Variation  with  Solar  Cycle 


The  figures  in  this  appendix  illustrate  how  the  diurnal/seasonal  pattern  of  occurrence 
statistics  (S4  >  0.5)  varies  as  a  function  of  SSN  at  Huancayo  (VHP),  Manila  (VHP),  and 
Ascension  Island  (L-band).  There  are  six  plots  for  each  station,  one  for  each  of  six  ranges 
of  SSN:  All  SSN,  SSN  <  50,  50  <  SSN  <  100,  100  <  SSN  <  150,  150  <  SSN  <  200,  and 
SSN  >  200.  In  all  plots,  the  data  have  been  limited  to  those  cases  for  which  KpSS  <  2°. 
The  bin  sisje  in  all  plots  is  30  minutes  in  GMT  by  30.5  days  (one  month)  in  day  of  the  year. 

I'h^i  heavy  solid  line  on  each  plot  shows  the  variation  of  the  time  at  which  both  El-region 
ends  of  the  IPP  field  line  with  day  of  the  year,  and  the  heavy  dashed  line  shows  the  valuation 
of  the  time  at  which  P2-iayor  sunset  occurs  (defined  as  x  =  107®).  The  horizontal  dashed 
lines  indicate  the  equinoxes. 
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Appendix  B.  Data  Plots  —  Variation  with  Geomagnetic  Activity  (Kp) 


The  figures  in  this  appendix  illustrate  how  the  diurnal/seasonal  pattern  of  occurrence 
statistics  (S^  >  0.5)  varies  as  a  function  of  geomagnetic  activity  level  (as  specified  by  Kp  at 
the  time  of  observation)  at  Huancayo  (VHP),  Manila  (VHP),  and  Ascension  Island  (L-band). 
There  are  four  plots  for  each  station,  one  for  each  of  four  ranges  of  Kp-.  All  A'p,  Kp  <  2o, 
2o  <  SSN  <  4o,  and  Kp  >  4o.  The  bin  size  in  all  plots  is  30  minutes  in  GMT  by  30.5  days 
(one  month)  in  day  of  the  year. 

The  heavy  solid  line  on  each  plot  shows  the  variation  of  the  time  at  which  both  E-region 
ends  of  the  IPP  field  line  with  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation 
of  the  time  at  which  F2-layer  sunset  occurs  (defined  as  ^  =  107*).  The  horizontal  dashed 
lines  indicate  the  equinoxes. 
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Appendix  C.  Data  Plots  —  Variation  with  Geomagnetic  Activity  (Kpss) 


The  figures  in  this  appendix  illustrate  how  the  diurnal/seasonal  pattern  of  occurrence 
statistics  (S4  >  0.5)  varies  as  a  function  of  geomagnetic  activity  level  (as  specified  by  Kp 
at  the  time  of  local  sunset)  at  Huancayo  (VHF),  Manila  (VHF),  and  Ascension  Island  (L- 
band).  There  are  four  plots  for  each  station,  one  for  each  of  four  ranges  of  Kpss-  All  Kpss, 
^  2o,  2o  <  SSN  <  40,  and  Kpss  >  4o.  The  bin  size  in  all  plots  is  30  minutes  in  GMT 
by  30.5  days  (one  month)  in  day  of  the  year. 

The  heavy  solid  line  on  each  plot  shows  the  variation  of  the  time  at  which  both  E-region 
ends  of  the  IPP  field  line  with  day  of  the  year,  and  the  heavy  dashed  line  shows  the  variation 
of  the  time  at  which  F2-layer  sunset  occurs  (defined  as  x  =  107").  The  horizontal  dashed 
lines  indicate  the  equinoxes. 
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